Published in The Astrophysical Journal, 539:L137-L141, 2000 August 20 
Preprint typeset using L^T^ style emulateapj v. 21/08/00 



EXTENDED [C l] AND "CO(5^4) EMISSION IN M17SW 

J. E. Howe', M. L. N. Ashby^, E. A. Bergin^, G. Chin\ N. R. Erickson', P. F. Goldsmith'*, M. 
Harwit^, D. J. Hollenbach'', M. J. Kaufman', S. C. Kleiner', D. G. Koch*, D. A. Neufeld', B. 
M. Patten-, R. Plume-, R. Schieder'", R. L. Snell', J. R. STAUFFER^ V. ToLLS^ Z. WANG^ G. 
Winnewisser'", Y. F. Zhang^, and G. J. Melnick^ 

Received 2000 May 1; accepted 2000 June 23; published 2000 August 16 

ABSTRACT 

We mapped a 13 x 22 pc region in emission from 492 GHz [C l] and, for the first time, 551 GHz '■'CO(5^4) 
in the giant molecular cloud M17SW. The morphologies of the [C l] and '^CO emission are strikingly similar. 
The extent and intensity of the [C l] and '^CO(5-^4) emission is explained as arising from photodissociation 
regions on the surfaces of embedded molecular clumps. Modeling of the '^CO(5^4) emission in comparison to 
'^CO(l^O) indicates a temperature gradient across the cloud, peaking to at least 63 K near the M17 ionization 
front and decreasing to at least 20 K at the western edge of the cloud. We see no correlation between gas density 
and column density. The beam-averaged column density of C I in the core is 1 x 10'^ cm"^, and the mean column 
density ratio A^(C l)/A^(CO) is about 0.4. The variations of N{C l)W(CO) with position in M17SW indicate a 
similar clump size distribution throughout the cloud. 

Subject headings: ISM: individual (M17SW) — ISM: molecules — ISM: structure — radio lines: ISM — 
submillimeter 



1. INTRODUCTION 

M17SW is a giant molecular cloud at a distance of 2.2 kpc 
• (Chini et al. 1980) illuminated by the M17 OB association 
' about 1 pc to the east. Studies of molecular excitation (Snell 
et al. 1984), molecular line profiles (Martin, Sanders, & Hills 

■ 1984), and maps of molecular emission (Stutzki & Giisten 
' 1990) from the M17SW core have all concluded that the dense 

gas is highly clumped. Early photochemical models of molec- 
' ular cloud/H II region interfaces predicted an overlying layer 
\ of ionized carbon and a thin underlying layer of neutral carbon 
within a photodissociation region (PDR) between the ionized 
and molecular gas components (e.g., Langer 1976a; Tielens & 
Hollenbach 1985). Given the edge-on geometry of M17SW 
and the H II region, observations of [C l] and [C ll] emis- 
sion along cuts through the interface region unexpectedly de- 
. tected these species more than a parsec into the molecular cloud 

■ (Keene et al. 1985; Genzel et al. 1988; Stutzki et al. 1988). 
Genzel et al. and Stutzki et al. concluded that far-ultraviolet 
(FUV) radiation penetrates the interclump medium and excites 
PDRs on the surfaces of illuminated clumps. Subsequent multi- 
component PDR models of the cloud have accounted for the in- 
tensity of the emission from numerous far-infrared and submil- 
limeter atomic and molecular lines (Burton, Hollenbach, & Tie- 
lens 1990; Meixner et al. 1992). Sekimoto et al. (1999) present 
observations of [C l] emission extending more than a degree 
southwest of the M17SW cloud core. 

The dual radiometers on board the Submillimeter Wave As- 



tronomy Satellite iSWAS) provide the ability to simultaneously 
observe the emission from the 7 = 5^4 rotational transition of 
'^CO at 550.926 GHz and the ^Pi ^ ^Pq fine-structure tran- 
sition of C I at 492.161 GHz. The roughly 4' beamsize of 
SWAS affords us the opportunity to map large regions in the 
'"'CO(5— >4) and [C l] lines and study the large-scale physical 
conditions in and structure of the dense interstellar medium. 
This is especially important for the '^CO(5^4) line, since few 
maps are available in the mid-/ '^CO transitions beyond the 
small-scale strip maps of Graf et al. (1993). Whereas the criti- 
cal density (wcr 2 x lO^* cm"^) and upper-level energy (£„ = 24 
K) of [C l] enable it to be easily excited wherever neutral carbon 
is present in the dense ISM, the '^^CO(5^4) line {n„ --2x10^ 
cm"-'; Eu = 79 K) will preferentially probe warmer, denser re- 
gions. 

We obtained with SWAS the first large-scale map of the 
'^CO(5->4) fine in M17SW, as well as the [C l] line over the 
same region. We combine these results with millimeter '^CO 
observations to probe the physical conditions in the extended 
cloud and the spatial and density structure of the bulk of the 
molecular gas. 

2. OBSERVATIONS AND RESULTS 

We mapped the [C l] and '^CO(5^4) fines in M17SW with 
SWAS between 1999 June 30 and 1999 July 12. We observed a 
21' X 35' region (13 x 22 pc) in mapping mode with a 1.'6 grid 
spacing, nodding the telescope between the source position and 
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Fig. 1. — Integrated intensity maps of emission from 492 GHz [C I] 
(left panel) and 551 GHz '^00(5^4) (center panel) in M17SW observed by 
SWAS, compared to a ''CO(l— »0) map smoothed to equivalent spatial resolu- 
tion (right panel; Wilson et al. 1999). Mapping offsets are relative to position 
a = 18''20'°22?1, 5 = -16° 12'37" (J2000). Contour levels for all maps range 
from 0.1 to 0.9 of the peak in steps of 0.1. Peak intensities (7^* scale) are 57.2 
K km s-i for [C I], 56.9 K km s"' for "CO(5^4), and 38.2 K km s"' for 
"CO(l^O), integrated over the Vlsr range 14-26 km s~'. 

a reference position at a = 18''25'"13?1, 5 = 
-16°50'37"(J2000). Integration times are typically about 3 
minutes, except for 16 central positions observed for H2O de- 
tections and described in Snell et al. (2000). The spatial resolu- 
tion of the observations is about 4', and the spectral resolution 
is about 1 km s"\ sampled with a spacing of 0.6 km s"'. The 
[C I] and '^CO(5^4) linewidths range from about 2-9 km s"^ 
(typically 5 km s"'), so the lines are well resolved. We reduced 
the spectra using the standard SWAS pipeline procedure, and 
removed small first-order polynomial baseline offsets. Details 
of the SWAS receivers, beam shapes, backend spectrometer, 
and observing modes are presented in Melnick et al. (2000). 
We present the maps herein on the 7^* scale, uncorrected for the 
SWAS main beam efficiency r?mb = 0.90. 

We show integrated intensity maps of the [C i] and 
'^CO(5^4) emission in Figure 1, as well as a '^CO(l^O) map 
from Wilson, Howe, & Balogh (1999), observed at 47" reso- 
lution with the 14 m telescope of the Five College Radio As- 
tronomy Observatory (FCRAO) but smoothed to 4.'0 resolution 
for comparison to the SWAS observations. The maps include 
emission within the velocity range (Vlsr) 14-26 km s~\ which 
excludes line-of-sight velocity components not associated with 
the Ml 7 H 11 region. The similarity of the morphology and ex- 
tent of the [C I] and '^CO(l— >0) emission is particularly strik- 
ing. This behavior of the [C i] and low-/ "CO emission has 
also been observed in other giant molecular clouds, notably 
Orion A (Tauber et al. 1995; Ikeda et al. 1999; Plume et al. 
2000), S140 (Plume, Jaffe, & Keene 1994), W3, L1630, and 
Cep A (Plume et al. 1999). The morphology of the "CO(5^4) 
emission is also similar to the [C I] emission but not as ex- 
tended, except at the eastern edge nearest the H II region. The 
"CO(5— >4) emission peaks slightly eastward of the [C I] and 
'^CO(l— >0) emission peaks, indicating higher temperatures or 
gas densities at the eastern edge of the cloud. We explore this 
further in §3. 

The shapes of the [C l] and '^CO lines are quite similar, with 
the [C l] and '^CO(l^O) line profiles being nearly identical. 
There are some subtle differences between the ^^CO(5^4) line 
and the other Unes, most notably a slight velocity shift at the 
map center position and slight variances in the red side of the 
line profile in the northeastern region. Overall, however, the 
spectra reveal a dynamical as well as morphological similarity 
between the [C i] and "CO gas components. 
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Fig. 2. — Model results for the distribution of '^CO colvmm density (left 
panel) and kinetic temperature (right panel) in M17SW. Position offsets are 
relative to o = 18h20'"22M, 5 = -16°12'37" (J2000). Contour levels range 
from 0.1 to 0.9 of the peak in steps of 0.1 forWC'^CO) (peak value 6.1 X lO"^ 
cm~^), and 0.3 to 0.9 of the peak in steps of 0.1 for Tk (peak value 63 K). 
Regions with "C0(5— >4) detections less than la are blanked. 

3. DISCUSSION 

3.1. Temperature, Density, and Column Density Structure 

We used the "CO(5^4) and "CO(l^O) integrated intensi- 
ties together with an iterative algorithm incorporating the Large 
Velocity Gradient (LVG) radiative transfer approximation (e.g., 
Scoville & Solomon 1974) to model the "CO column density 
N(}^CO), H2 number density ^(Hi), and kinetic temperature Tk. 
In the LVG modeling we used colhsion rates for CO with para- 
H2 (Flower & Launay 1985). Substituting ortho-Hi rates re- 
sulted in insignificant differences in the model results. For each 
position, we began with initial estimates of Tk= 10 K, A^("C0) 
appropriate for LTE, and n(H2) from Snell et al. (1984) for the 
core or n(H2) = 1 x 10^ cm~^ elsewhere (Wilson et al. 1999). 
We used these estimates as input for the LVG algorithm, to- 
gether with velocity widths derived from the "CO(l^O) data, 
to produce model "CO(5^4) and "CO(l^O) integrated in- 
tensities. We iterated the model with density as the free parame- 
ter until the model "CO(5^4)/"CO(l ^0) integrated intensity 
ratio reproduced the observations. If the 7 < 5 level populations 
thermalized with too low a ratio (i.e. further increases in density 
no longer affected the line ratio), we incremented Tk and reit- 
erated until the observed ratio was reached. We then iterated 
the model with A^("CO) as the free parameter until the abso- 
lute intensities of the lines agreed with observations. This value 
of Af("CO) was then used again in the iteration with n(H2) to 
check that Une opacity effects did not change the line ratio, and 
the entire loop repeated as necessary until the model Une inten- 
sities and ratios agreed with the observations. This procedure 
then yields the minimum temperatures necessary to explain the 
fine ratios. 

The model results show clear evidence that the molecular 
cloud is heated by the H II region/ionization front at its eastern 
edge. In Figure 2 we show the temperature distribution over 
the region mapped in "CO(l^O) as well as the derived "CO 
column density distribution. The peak beam-averaged column 
density is 6.1 x 10^^ cm"^, with a mean value of 1.2 x 10^^ 
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Fig. 3. — (Top) The linear correlation of the integrated intensity (Tyss scale) 
of [C I] with "CO(l^O) at each position in M17SW. The dashed line traces 
the least-squares linear fit to the data (slope 0.84). The dotted line plots the 
equivalent fit for the Orion A cloud (slope 0.5; Ikeda et al. 1999). (Bot- 
tom) The derived ratio of the column densities of C I to CO [where Af(CO) 
= 50A'(''CO)] is plotted against the "CO column density at each position 
in M17SW. Errorhars are ± Ict as described in the text. The dashed line 
marks the mean value of N(C l)/N(CO) (0.37) for all data points. The iii- 
angles plot the value of ^(C I)/7V(CO) in bins containing equal numbers of 
'^CO molecules, with the abscissa at the center of each binning range. 

cm"^. The column density map appears virtually identical to 
the ^^CO(l^O) integrated intensity map (Fig. 1), a property 
attributable to the low opacity of the ^^CO lines (r < 0.6 for 
either line). The highest gas temperatures are found along the 
eastern edge of the cloud, peaking at about 60 K near the cloud 
core. The mean cloud temperature is 35 K, similar to the global 
cloud temperature of 30 K derived by Wilson et al. (1999). The 
temperature we derive for the core region is also consistent with 
the temperatures derived from millimeter observations of opti- 
cally thin CH3CCH emission (Bergin et al. 1994). Our obser- 
vations apparently are insensitive to the high temperature gas 
(Tk > 200 K) responsible for CO(7— »6) line emission observed 
in the H 11/H2 interface region, however this gas is estimated to 
comprise only about 5 percent of the total column density (Har- 
ris et al. 1987) and would contribute only about 10 percent of 
the observed peak ^^CO(5— >4) intensity. 



We derive molecular gas densities in the range 1-6 xlO^ 
cm"^ over the entire region mapped, with a median value of 

3 X 10^ cm~^, about half as large as densities derived from 
multi-transition CS observations of the core region of M17SW 
(Snell et al. 1984; Wang et al. 1993). In contrast to the well- 
ordered column density and temperature distribution across the 
M17SW cloud, the density distribution appears almost random. 
In particular, there is no correlation between the gas density 
and the column density. Instead, high- and low-density regions 
are peppered randomly throughout the region. Since the excita- 
tion of the '^CO(5^4) emission depends both on the tempera- 
ture and density, with higher densities requiring lower temper- 
atures for the same 7 = 5 fractional population, we investigated 
whether there was any trend between Tk and n(H2) as an arti- 
fact of the modeling algorithm. We see no correlation between 
the gas temperature distribution and the density distribution in 
the model results. 

Although the model distributions of temperature and density 
are uncorrelated, the densities we derive are sensitive to the 
model cloud temperature. For example, uniform cloud tem- 
peratures of 50 K and 100 K yield median cloud densities of 

4 X lO'' and 1 x 10* cm~^, respectively. The model column den- 
sities are insensitive to temperature, decreasing by only 20 per- 
cent if Tk = 100 K. Since such large column densities of hot 
gas are difficult to produce throughout the cloud, and since our 
model mean cloud temperature is already somewhat higher than 
the Wilson et al. (1999) global cloud temperature, we consider 
4 x lO'* cm~^ as a reasonable lower limit to the mean cloud den- 
sity. 

The modeling results for the density and coluimi density of 
the molecular gas in M17SW require that the gas is highly 
clumped. We derive the volume filling factor of the gas by di- 
viding the ratio of the peak H2 column density (4 x 10^^ cm"^, 
assuming a '^CO abundance of 1.5 x 10"^ relative to H2) to 
the mean H2 number density (3 x 10^ cm~^) by the Unear ex- 
tent of the core (3 pc core size, deconvolved from the SWAS 
beam). From these quantities, we derive a volume filling factor 
of 0.02. Stutzki & Gusten (1990) derive toward the central 13" 
of the core a volume filling factor of 0.13 for gas with densi- 
ties greater than 10^ cm"^, reflecting a higher concentration of 
clumps at the core peak region. 

3.2. The Correlation o/[C i] and ^^CO Emission 

The large extent of the [C i] emission in M17SW and its sim- 
ilarity to isotopic CO emission was first noted by Keene et al. 
(1985) and Genzel et al. (1988) along linear cuts through the 
ionization front/molecular cloud interface. The large-scale [C l] 
and '^CO maps in Figure 1 show that this similarity extends 
over the entire region mapped. We plot in the top panel of Fig- 
ure 3 the observed correlation of the integrated intensities of the 
[C l] and '^CO(l^O) emission. We converted the data to the 
main-beam temperature scale T^b to remove differences in tele- 
scope beam efficiencies from the observed relation {r]mb = 0-52 
for the FCRAO data; Wilson et al. 1999). The observed trend is 
well fit (linear least-squares, correlation coefficient 0.95) by the 
relation / WV([C l]) = (0.84 ±0.02) / TrnbdVC^CO(l^0))~ 
(0.4 ±0.4). The slope of this relation is similar to the slope ob- 
served between [C l] and '^CO(2^1) in the Orion Bar region 
by Tauber et al. (1995; ^ 0.8), but is somewhat higher than the 
relation reported by Ikeda et al. (1999; ~ 0.5) for the Orion A 
cloud as a whole, and for the ensemble of clouds W3, L1630, 
S140, and Cep A (Plume et al. 1999). We propose a possible 



4 



[C I] AND i3cO(5^4) EMISSION IN M17SW 



explanation for this in §3.3. 

We used the SWAS [C i] observations to model the C i col- 
umn density A^(C i) in a manner similar to the '^CO LVG 
modeling, using the collision rates of Schroder et al. (1991). 
Since the spatial distribution and velocity profiles of the C I and 
^^CO hues are nearly identical, we used the temperatures and 
densities derived from the ^^CO observations, along with the 
'^CO(l— >0) linewidths, as fixed input parameters to the LVG 
algorithm. We then iterated the models with N(C l) as the sole 
free parameter until the model |C l| intensities matched the ob- 
servations. NiC I) ranges from 1.1 x 10'* cm"^ to 1.1 x 10'^ 
cm~^ with a median value of 1.7 x lO'^ cm~^. The [C i] line 
is optically thin, with a maximum opacity less than 0.6 at the 
core. Our peak value of A^(C l) is about half that obtained by 
Sekimoto et al. (1999) in a 2' beam. The spatial distribution of 
A^(C I) is nearly identical to the distribution of A^('^CO) shown 
in Figure 2. 

The derived C I column densities show an interesting behav- 
ior when compared to the '^CO column densities. In the bottom 
panel of Figure 3 we plot the ratio of A^(C I) to the CO colunm 
density N{CO) (calculated assuming a CO/'^CO abundance ra- 
tio of 50) as a function of '^CO column density for each model 
position. Errorbars reflect the random errors of the [C l] and 
'^CO(l— >0) intensities, but not systematic uncertainties in the 
LVG models. Whereas at the low column densities at the edges 
of the molecular cloud the A'^(C i)/NiCO) ratio ranges from 0.1 
to 1.0, at high column densities in the core of the cloud the 
lower and upper bound for the range asymptotically approach 
the mean N(C l)/Af(CO) ratio for all positions (0.37). This is 
similar to the behavior exhibited for the Orion A cloud, except 
the mean column density ratio is much lower, only about 0. 10 in 
the Orion core regions (Plume et al. 2000). The behavior of the 
N{C l)/N(CO) ratio is difficult to understand in the context of 
chemical models where the neutral carbon is produced in situ 
in the molecular gas (e.g., Langer 1976b; Pineau des Forets, 
Roueff, & Flower 1992), and would thereby predict a more or 
less constant column density ratio. Likewise, simple plane par- 
allel FDR models (e.g., Langer 1976a; Tielens & Hollenbach 
1985; Hollenbach, Takahashi, & Tielens 1991) predict a rising 
value of A^(C i)/A^(CO) with decreasing molecular gas column 
density, but not the increased scatter we observe in the ratios. 
Remarkably, summing the A^('^CO)-ordered data of Figure 3 
into bins of equal total number of '■'CO molecules (in this case, 
2 X 10'^ cm"^ X Abeam^ whcrc Abeam is the area subtended by 
the SWAS beam at the distance of M17SW), and then comput- 
ing A^(C l)/A'(CO) for the binned data, shows that A'(C l)/A^(CO) 
is independent of the mean '^CO column density for each bin 
(see Fig. 3). The resolution to the behavior of N{C i)/N{CO) 
lies in the structure of the M17SW molecular cloud. 

3.3. PDR Emission and the Structure ofMlJSW 

Theoretical models of PDRs are quite successful at predict- 
ing the observed '^CO(5— >4) and [C I] line intensities and the 
column density of C ifor the physical conditions in the region 
mapped by SWAS. Based on the arguments of Stutzki et al. 
(1988) and Meixner et al. (1992), we expect the incident FUV 
field within the cloud to range from ^ 10 Go at the northwest- 
ern edge to 4 X lO'* Go at the ionization front located about 1 .'5 
northeast of the map reference position, where Gq is in units 
of the Habing flux. Our LVG modeUng indicates a fairly nar- 
row range of densities, 1-6 x 10^ cm~^. PDR models predict 



a column density A^(C l) ~ 4 x 10'^ cm"^ insensitive to varia- 
tions in Go and only weakly dependent on density (Tielens & 
Hollenbach 1985). This is similar to the mean value of A^(C I) 
derived for the M17SW cloud. For the peak A^(C i) we derive, 
only a few PDRs along the line of site through the core are re- 
quired. Kaufman et al. (1999) present updated homogeneous 
PDR models covering the range of densities and FUV fluxes 
appropriate for M17SW, and have recently incorporated '^CO 
into the models for comparison with SWAS observations. These 
models successfully predict the same range of '•'CO(5^4) in- 
tegrated intensities, for the densities and FUV fluxes we find 
in M17SW, as those observed by SWAS. The PDR models of 
Koster et al. (1994) also predict the peak '^CO(5^4) fine tem- 
perature we observe, for Go = 10"* and n(H2) = 10^ cm"^. The 
Kaufman et al. (1999) models match the observed [C l] inte- 
grated intensities for regions of the cloud with Go < 100, but are 
a factor of about 3-4 too low for the core of the cloud. Again, 
this is easily explained by requiring several PDRs along the line 
of sight through the core. 

The extent of the [C i] emission and the relationship of A^(C i) 
with the molecular gas coluimi density require not only that the 
gas is clumpy, but that the distribution of clump sizes is similar 
throughout the cloud, from the core out to the periphery. In a 
clumpy cloud penetrated by FUV radiation, PDRs arise on the 
surfaces of the dense clumps, with a thin shell of C i encom- 
passing the CO in the clump. Consider unresolved clumps of 
similar density, as is the case for the SWAS observations (see 
Stutzki & Glisten 1990). The beam-averaged column density 
A^(C I) depends on the clump surface area while N(CO) de- 
pends on the volume, so the ratio A'^(C i)W(CO) varies as the 
inverse size of the clumps. Thus, one might expect for a cloud 
with a range of clump sizes a large scatter in the column density 
ratio, as we see in the low column density regions where only 
a few clumps are probably subtended by the SWAS beam. As 
the beam samples larger column densities, i.e. more clumps in 
the beam, the clump size distribution is more completely sam- 
pled and the scatter in A^(C i)/N(CO) decreases, approaching 
the value appropriate for the mean surface area to volume ratio 
of the clump size distribution. By summing up the low col- 
umn density regions we increase the clump sample size, and 
we see from Figure 3 that the A^(C l)/N(CO) ratio for equal to- 
tal gas masses gives the same value of NiC i)/NiCO). Thus, 
the size distribution of clumps must be the same throughout 
the cloud. The diameter of a spherical clump with a C i/CO 
number ratio of 0.37 and H2 density 3x10^ cm"' is 0.04 pc, 
where we assume all carbon is either in C I or in CO and 
calculate the thickness of the C ilayer from Hollenbach et al. 
(1991). This diameter scales inversely with the H2 density. 
The lower A^(C l)/N(CO) ratio found in Orion is probably a 
consequence of a larger mean clump size, since A^(C l) is in- 
sensitive to changes in density and FUV flux. Since the inten- 
sity of optically thin fines scales directly with column density, 
observationally a lower A'^(C i)/NiCO) ratio produces a lower 
slope in the linear correlation of the [C l] fine intensity to the 
'^CO(l^O) intensity (e.g. Fig. 3). 
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grants 50 0090 090 and 50 0099 Oil. 
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